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The main application of Heavy Quark Effective Theory (HQET) and of Soft Collinear Effective 
Theory (SCET) is in establishing factorization theorems for exclusive and semi-inclusive decays of 
heavy mesons. However, the calculation of the soft factors from the HQET or SCET factorization 
relations is, as usual, impeded by the non-perturbative aspect of the strong interactions. In the 
hope of gaining some insights into some of these quantities we compute them in the 't Hooft model. 
We find that the B-meson shape function is exactly given by the square of the _B-meson light-cone 
wave-function. The structure of the B — tt structure function is more complicated: it is given by the 
product of wave-functions or by a resonant sum depending on the kinematics. The result simplifies 
dramatically in the chiral limit, where it can be compared with general arguments based on Heavy 
Meson-Chiral Perturbation theory. No attempt is made to use these results for applications to 
. phenomenology. 
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I. INTRODUCTION 

One of the most celebrated applications of Heavy Quark Effective Theory (HQET) 0,0, 0010] is to the calculation 
of the rate for the semileptonic decay of a B-meson. It is found that the differential decay rate (with respect to a single 
kinematic variable) is systematically predicted as an expansion in inverse powers of the heavy quark mass, l/m&, with 
Q\ \ the leading term completely fixed by symmetry, and hence calculable 0, @- In the case of charmless semileptonic 
I/"") • B-meson decays the systematic expansion breaks down in the "end-point region," namely, when the leptonic energy 
is as large as it can be. For x = 2E /m,B ~ 1 — (A/to^) 2 convergence is hopeless, while for 1 — x of order unity the 
expansion is very good. There is a region, 1 — x ~s A/ tub, however, for which all orders in the HQET expansion are 
equally large. The most singular terms in the expansion can be formally re-summed into a so-called "shape-function" 

Since the shape function involves the sum of an infinite number of terms, each characterized by an unknown constant 
\ (the matrix element of a local operator between B- meson states), it is largely unknown. Some basic properties can 
be readily derived {e.g., normalization follows from current conservation), but in practice this is far from sufficient 
to determine it with the precision required for applications. The favored solution to this impasse is to show that the 
shape-function is approximately universal and then (approximately) cancel it from measurements of, say, charmless 
semileptonic and radiative -B-decays. 

In this short note we propose to compute the shape function in a the large iV c -limit of QCD in 1 + 1 dimensions, 
the 't Hooft model. While the model can not be used to compute parameters of phenomenological interest, it has 
been extensively studied in the past and has given remarkable insights into the dynamics of the strong interactions. 
■ In the context of heavy mesons, calculations in the 't Hooft model gave some of the earliest indications that the heavy 
quark expansion converges rather more quickly for some quantities than for others: relations among form factors for 
B — > Dlv and the predicted normalization of these form factors at zero recoil hold more accurately than relations 
among heavy meson decay constants Il2| . Some results are interesting but are unlikely to have a counterpart in 
four dimensions: in Ref. [Hj it was shown that in the chiral limit the form factors for the decay B — > -kIv are exactly 
given by a single pole formula (and the residue is fixed by symmetry considerations). And sometimes the results have 
served to raise warning flags about possibly unjustified assumptions in the four-dimensional analysis. For example, in 
contrast to what has been argued informally 14|, a numerical solution to the model 15] shows that the l/mt, expansion 
of the lifetime of the i?-meson contains corrections of first order in the expansion parameter, that is, (l/mf,) n with 
n = 1. To be sure, this result is controversial. It is trivial to show analytically that the (l/mb) n , n = 1 corrections are 
absent in the chiral limit 0], but attempts to extend this result to the non-chiral limit (where the result of Ref. 0] 
applies) are far from rigorous. Moreover, the numerical result in is consistent with the theoretical observations 
in Ref. |l7l Il8| that the l/m& expansion for the smeared B-width has no (1/mf,) 1 term. In sum, the 't Hooft model 
is a suitable tool for testing proposed methods and mechanisms in 4-dimensional QCD that should just as well apply 
in 2-dimensional QCD, but it is not likely a good model for phenomenological applications. 
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There is another quantity of interest we will compute in this note. The shape-function can be described as the 
expectation value in a £?-meson state of a non-local b-quark bilinear operator, with a Wilson line between quarks to 
ensure gauge independence. A similar quantity of interest is the matrix element between different states, say, a B- 
meson and a w- meson, of a non-local quark bilinear (with properly chosen flavor quantum numbers). This generalized 
parton distribution (GPD) plays a central role in the rece ntly establishe dll9l so ft pion theorems for certain soft factors 
in the Soft Collinear Effective Theory (SCET) [U H H Hi H |U |H EfiM, El SEHSl description of exclusive 
B decays to final states containing, possibly, several soft pions. This, in turn, is useful in extending the app licability 
of SCET to the calculation of non-resonant final states, as in, for example, B — » Ktt££ [33] or B — > Ktt^\34. l35l| - 

The paper is organized as follows. We begin in section [H] with a review of the 't Hooft model. This is intended to 
establish notation and to review some techniques, but it is not recommended as a place to learn about the model, 
which can be best done by going to the original sources [3(1 l37l EM 13^. Section ITTT1 shows that the 't Hooft wave- 
fucntions for a meson are, up to normalization, the light-cone wave-functions for that meson. This is an old result, 
but we use it here as as simple example of the calculational techniques used in Sees. HVl and IVl where the shape and 
soft functions are computed, respectively. In Sec. we endeavor to determine the chiral limit of the shape functions 
obtained for general values of the quark masses. This requires some knowledge of Heavy-Light form factors, and of 
how the chiral limit is approached for these form factors, so we have included a review of those results as the first 
sub-section in Sec. [V] and we have separated the chiral limit of the soft functions into a separate sub-section as well. 

Section IVII is the most interesting: we discuss the interpretation of the results just obtained and compare them 
to the literature. We find that the method for computing soft functions proposed in Ref. 0] is not quite correct. 
Moreover, we propose a simple (but incomplete) fix. And we find that the shape function is simply the square of the 
light-cone wave-function. These two resulsts are tightly connected in the 't Hooft model and, as we remark in the 
concluding section fVIII it would be interesting to determine if this connection persists in 4-dimensional QCD. This 
would lead to interesting phenomenology. Finally, some technical issues, concerning scaling functions, have been left 
to an appendix. 



II. REVIEW OF THE MODEL 

This model has been extensively studied and our work relies on technology pioneered by 't Hooft[3(|, Callan, Coote, 
and Gross 37], and Einhorn 38J. In these papers the bound state equations were derived; and it was shown that the 
scattering amplitudes — and the form factor in particular — can be written entirely in terms of interactions among the 
meson bound states, with no quarks in the spectrum or in the singularity structure of the amplitudes. 

We recall the features of the model which make it solvable, and refer the reader to the original papers for details. 
The dynamics arc defined by the Lagrangian, 

C = —TtF^F^ + X^„(V(i0 M - goAJ - m a )^ a , (1) 

a 

where is an SU(N C ) gauge field, is its field strength and ip a is a Dirac fermion of mass m a . In the large- N c 
limit, the gauge coupling is scaled with N c : g 2 — g^N is held fixed as N c — > oo. The label a runs over flavors of quark, 
with bare masses m a . 

The theory is most conveniently quantized in light-cone gauge. Because there are no transverse dimensions, set- 
ting A_ = eliminates the gluon self-coupling. It also serves to project gamma matrices onto a single component 
in any Feynman graph that has just gluon vertices and (— ) component current insertions on fermion lines. In this 
gauge the gluon propagator is the inverse of d_. The infrared divergence in the gluon propagator, i/k_, is regulated 
by taking the principal value at the pole. 

Spinors arc conveniently split into L and R components, 

ip R = P R ip = !7+7_-0, = P L t/j = ±7_7 + ?/>, (2) 

and in A_ = gauge i^l is not independent, 

^L = - j f^^R- (3) 

The leading term of the 1/N C expansion is the sum of planar graphs. The equations for the full propagator and 
self-energy can be solved exactly, with an extremely simple result: the net effect of all gluons starting and ending on 
the same fermion line is just to renormalize the quark mass appearing in the propagator, 



m l^rn 2 a = m 2 a - g 2 /n, 



(4) 
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FIG. 1: Feynman diagram representation of the full meson-quark vertex (fc)r^ 12 ^ (p, k)S^ 2 ' (k — p) (propagators have been 
included so the vertex function is amputated). The cross denotes an insertion of the current rpi'y-fa. 



so the full quark propagator is 



S(«0(k) = -. (5) 

w k 2 - m 2 + ie w 



After making this shift, all remaining gluon interactions enter as ladder-type exchanges. Crossings would require 
either gluon self-couplings, which are absent, or non-planar graphs, which are higher order in l/N c . The Yang-Mills 
coupling constant g has dimensions of mass, and we choose units such that g 2 /ir — 1, leaving m 2 as the dimensionless 
numbers parametrizing the theory. As is well known, there is a discrete spectrum of free mesons. 

The full meson-quark vertex, T(p, k) (see Fig.^J, satisfies a Bethe-Salpeter-like equation: a gluon exchange between 
incoming and outgoing quarks leads to 



r< 12 ) (p, k) - 1 = -- 4 d\ 7 — V^fS (1) (?)r (12) (p. i)s (2) (q-P). (6) 



7t J (q- — k 

Here and below, the mark through the integral sign reminds us to take the principal value of the integrand. Eq. © 
shows that T(p, k) only depends on k through the variable x = k-/p~, so we define T(p 2 , x) = T(p, k). The equation 
for the vertex function is solved in terms of the solutions of the associated eigenvalue problem, which has the physical 
interpretation of the bound state equation: 

^ n (x) = ( M + ^L) Mx) - f j^^y). (7) 
\ x l-xj J {y — x) 1 

The function (j) n (x) is commonly referred to as the 't Hooft wave- function of the n th eigenstate, with mass /i„, 
and x = k-jp- is the fraction of light-cone momentum carried by the out-going quark. The out-going and in-going 
quarks may have different flavor, and hence their masses mi and mi are in general distinct. When needed we will 

(12) 

explicitly denote the wave-function by the flavor of the quarks, as in <jr n . It is conventional to choose a standard 
normalization, J Q dx 4> n (x)4> m (x) = 5 nrn . The <fi n (x) vanish at the boundaries, and consistency of (JJJ requires that as 



x 



0, 4> n {x) ~ x 131 , with 



7r/3i cot 7r/3i = 1 — m\ , (8) 



and similarly with the replacements m\ — > mi and x — ► 1 — x as x — >l,as dictated by the boundary behavior of 
the Hilbert transform. The bound state equation does not have solutions in terms of known functions, but may be 
readily solved numerically. 

The range of x for the bound states is always in the interval [0, 1], and <f) n — outside of this range; but Q 
determines as well the full meson-fcrmion-antifcrmion vertex, 

for values of z ^ [0, 1]. This includes z complex, corresponding to the general case where one or more of the lines of 
the meson-quark vertex is not on-shell in its physical region. $„(z) is analytic in the complex plane, with a cut on 
the real axis from to 1. When x € [0, 1], <fr n (x) is defined by the principal value prescription, and 

= f-/4 + — + Mx), (io) 

V x 1 — x I 



FIG. 2: Four quark scattering amplitude. 



in accordance with l(7|. Since <p n (x) is finite, <J> n (x) has zeros where the first factor on the right vanishes, and these 
are the values x± where the quarks would be on-shcll. These zeros of the vertex function cancel quark poles in the 
propagators of loop amplitudes to ensure that no quark singularities appear in gauge-invariant Green functions. 

All loop integrations are simplified by the fact that <fin(x) is a function of x = k-/p- only and is independent of 
p+. When wave functions and propagators appear in loop integrals, only the latter depend on p+, so the f dp+ is 
over rational functions and can be computed explicitly by contour integration, leaving a single integral over one real 
variable. 

The full meson-quark vertex is now easily expressed in terms of the eigenfunctions, 

iv,x) = i-]r^M, (ii) 

„ P 



where 



1 

dx<t> n {x) . (12) 







In the interval < x < 1 a handy alternative expression is available: 



2 / 2 _m?_^ E c^ 
The four quark connected scattering amplitude can be similarly obtained: 

2*L [ S x * 2 -/4 [ ' 

Here the kinematics is as indicated in Fig,[2]and we have defined s = p — k, x = r_/s_ and z = p_/s_. 

Below we will compare some of our results with those obtained from general considerations using heavy quark and 
chiral symmetries. It is therefore useful to present the 't Hooft wave-functions in the limiting case of one heavy quark 
(and a light anti-quark). In heavy quark effective theory (HQET) Green functions are expanded about the large 
momentum of the heavy quark. If the momentum of the heavy meson of mass fj, n is p, then v = p//i n is its velocity. 
We take mi to be the heavy quark mass, ra\ 3> 1 while keeping m 2 fixed, and write p = m\v + p and k = ni\v + fc; 
see Fig.^ Then the argument of the 't Hooft wave-function is x = k-/p- = 1 + — f>_)/mi«_ + • • • and the large 
mass limit is obtained by writing a function of fe_ = (1 — x)miv^: 

i>n(k-) = 0„(1 - k-/miV-) (15) 

(we have dropped the tilde since fc_ is a dummy variable). Note that fc_ — p_ = fc_ —p~ is the light quark momentum 
going into the graph. The restriction x = k-jp- < 1 is then the same as the statement that the wave-fucntion is 
non- vanishing only if the light quark momentum points out (when the heavy meson momentum points in). 

Making this change of variables, writing \x n = m\ + A n and taking the mi — > oo limit, the 't Hooft wave-function 
takes the form 



2A„^(fc_) = ( — + Mk-) - v- r , d9 ~ (16) 



It is easy to check that this is the equation one obtains starting directly from the HQET Lagrangian^ 
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III. LIGHT-CONE WAVE-FUNCTION 



Define the light cone wave- functions for a meson M n by 

h(0 =~T f ^ e- <(1_6)x ~M0|Wt(z-)|M n (p)) (17) 

In J 27T 

where the gauge invariant non-local operators W± , defined by 

W±(x) = ^{x)P[e^o A ] 7±75 ^(0) (18) 

have been evaluated along the x + = line. In the definition of W±, J A — J dy^ A fl (y), so in our computations, which 
are performed in light cone gauge, this is Jdy + A~ . The decay constant /„ is defined through 

(oi^ ( o)iM n(P )) = ( w» f * n I s ev ;;/ odd pari f!' d9) 

I i£p,vV ]n it n is odd (even parity) 

The antisymmetric two index tensor makes two dimensions special: it allows us to use an axial current as 
interpolating field for scalar mesons (and, similarly, to use a vector current for pseudo-scalar mesons). It is an easy 
exercise to compute the two-current correlator from which it follows that the decay constant is related to the integral 
of the wave-function in Eq. (|12fl , 

fn = \[^C n (20) 

The light-cone wave-function 0_ is given, up to normalization, by the 't Hooft wave-function. It is useful to review 
the derivation of this fact since the procedure is the same as in the more involved calculation of other soft functions 
below. Consider 

d 2 y e- ip - y (0\T[W(x-)^-ij(y)}\Q) . (21) 

In light-cone gauge this is really a three point function. To see this, take the path in the integral in Eq. (|18|) to be a 
straight line of y + = constant, from y~ = to y~ — x~ . After some spinor algebra, we have 

W-{x-)=tf R {x-)1> R (Q), 
W + (x~) = 4(ar)tfi(o) . 

Therefore, the correlator in Eq. (|21(l is related to the diagram in Fig.^by Fourier transformations. In order to perform 
the computation a simple trick proves useful. From the three point function 



T(p, k) = jd 2 x d 2 y e-^- +l ( fc -rt^(0|r[^(2/)Vfl(0)VU fl (x)]|0) , (23) 
which is a non-gauge independent quantity given by Eq. llljl . determine 

d 2 ye^>y(0\T[^y)^ R (0))\M n (p)) , (24) 



by extracting (and properly normalizing) the residue of the pole for the n-th meson state \M n (p)). At this point 
one can Fourier transform back to x space by integrating over J d 2 k e~ zk ' x and obtain (0\W(x~)\M n (p)} by setting 
x + = 0. This is the gauge invariant integrand of the definition of the light-cone wave- function in Eq. 117fl . The 
light-cone wave-function is the the Fourier transform of this with respect to x~ (which undoes one of the previous 
Fourier transforms). 

Let's see how this works, explicitly. Extracting the pole, expression l|23|) is 



d 2 y e« k ^(0\T[tl>tiv)MmMM = ^^<Mfl >2 ^ + • 7Z \, ? ~2+- ( 25 ) 

m{ + le (k — py — 77^2 + ie 
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where £ = fc_/p_. Following the strategy outlined above, we invert the Fourier transform and set y + = 0, 



wnrt <n*.«»H - / (g, «-"-'->»- i f g «-(e t p-^ +fc(t -^:^. +<€ < 26 > 

The integral over & + can be performed explicitly, by appropriately choosing a semicircle at infinity in the complex fc + 
plane to close the contour of integration, 

(0|T[^(y-)^(0)]|n) = f~ ^ e -*(*— _J (27) 

• /U Mm I — 

Since the argument of <!>„ is in the unit interval we can use Eq. (|10|l to relate it to the 't Hooft wave-function, 

/ AT /■! 

(0|T[4(y-)^(0)]|n) = -ip_y ^ jf d£ e* 1 ^-^^) (28) 

The coefficient in front is related to the decay constant by Eq. I)20|l . and the integral over £ is un-done by the Fourier 
transform that defines the light-cone wave function, Eq. (|17H . The result is zero for £ < or £ > 1 and, in the unit 
interval, 

^-(0=^(0/^- ( 29 ) 

As a fairly trivial check of the calculation, note that d£ </>„(£) = 1. As advertised, up to the normalization factor, 
the light-cone wave-function is the 't Hooft wave-function. This, of course, is a well known fact, but the method 
exemplified here with this calculation is precisely what we will use to calculate the shape and soft functions. 

The computation of the light-cone wave-function <j) + proceeds in an entirely analogous manner. It follows from 
Eq. J2J that the result is obtained by including in the computation above an additional factor of m\m2 / '4pi_j>2- > 
where pi are the quark momenta. Therefore, 

= - 2p 2™ U0/Cn ■ (30) 

Again we can check this by computing its integral. We need the useful relation 

"dx^J^™- = (-!)« £c n , (31) 
x(l — x) 



from which it follows that 

'0 

On the other hand, the integral of Eq. (|17|l gives 

1 i 



f 

Jo 



df = (-l) n+1 7^ = (-l)' i+1 ^ • (32) 



2p 2 _ v ' p 



d^„+(£) = r (0\W+(0)\Mn(p)) (33) 

P-Jn 



which agrees with Eq. 132fl after use of Eq. (|19|l . 



IV. SHAPE FUNCTIONS 

The shape functions are defined bv[^. IToj] 



/r n) (0 = ^ [dx-e-^-v-iBMbWPijK^TbmBnip)) _ (34) 



c+=0 



The notation suggests the meson, B, has a heavy &-quark and a light u-quark as constituents. In fact, our computation 
is valid for arbitrary masses of the quarks, but the notation was adopted since it is in the context of heavy mesons 
that shape functions often arise in practice. 
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FIG. 3: Non-resonant (left) and resonant (right) contributions to the Green's function for the soft function. The crosses denote 
insertions of the b R u R = V'Ik^.r current (or it's hermitian conjugate), both external lines are 6-quarks and the internal line is 
a it-quark. 



It is convenient to specialize to the case Y = 7_. Then, using \x = — vector-currents to interpolate for the mesons, 
we are led to consider a correlator of right-handed quark fields, 

Jd 2 xd 2 yd 2 z e *p(*-*)+**»<o|7 T [6t l (y)6K(0)6j i (ar)uj l (a;)uJ i (2r)6«(2r)]|0> . (35) 

As in the previous section, this quantity can be readily computed but is not gauge- invariant. However, undoing the 
Fourier transform over y and specializing to the line y + — does give a gauge independent quantity. 

Two distinct contributions to l|35l) arc shown in Fig. EJ It is easy to check that the second graph, involving the full 
scattering kernel, gives a vanishing contribution to the shape function. The first graph gives 

N C Y {12) (p, fc)r (2l) (-p, k - p)S (1) (fc)S (1) (k)S {2} (k-p). (36) 

We have indicated with superscript the quark flavors in these functions, with £=1,2 denoting q — b,u, respectively. 
For the vertex function the superscript (12) indicates an outgoing 6-quark and an incoming it-quark, and we adopt the 
same notation for other related quantities, e.g., the 't Hooft wave- functions. Applying LSZ reduction and integrating 
over k + this gives 

0l 15) (o4 2I) (i- omw-o, (37) 

where £ = fc_/p_. We have kept the subscript n arbitrary, indicating that the result is more general than we set out 
to obtain, namely, it is valid for the shape function of any of the states in the tower of which the B meson is the 

(21) (12) 

ground state (corresponding to n = 0). Using 4> n '(1 — £) = 0„ ; (£) we finally have 

/i n) (0 = (^ l5) (0) 2 (38) 

where it is understood the support is for < £ < 1 and the subscript in / reminds us of the choice Y = 7_ . 
From the normalization of the 't Hooft wave-functions it follows that 

[ de/i n) (o = i (39) 

Jo 

upon integration of Eq. (|38|l . This is a minimal test on our calculations since the normalization of the shape function 
follows from charge (6-number) conservation. 

Finally we express this result in terms of functions in the heavy quark limit. It is customary to write the shape 
function in terms of the residual momentum k- of the 6-quark. Our expression for the heavy quark 't Hooft wave- 
function, Eq. i|15f) . has the light quark momentum as argument. To write a heavy quark limit formula for the 
shape function as a function of the 6-quark residual momentum fe_ we use A n = ji n — mi, so that £ = fc_/p_ = 
(miu_ + k-)/[m\V- + A„t>_) = 1 + (fc_ — k n v-)/miV- + ■ ■ • . We obtain 

/i n) (M = [V>„(A„-fc_)] 2 (40) 

Note that this has support only for fc_ < A„. 
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FIG. 4: Non-resonant (left) and resonant (right) contribution to the Green's function for the soft function. The crosses denote 
insertions of the b R d R = ip\ R tp 3R current on the left and the d R u R = ip\ R ip2R current on the right. The outgoing external 
fermion is a 6-quark, the incoming external line is a u-quark and the internal line is a d-quark. 



V. "HEAVY-LIGHT" SOFT FUNCTION 



A. "Heavy-light" Form Factors 



In preparation for our study of the soft function we review some basic results for the transition current form factors. 
We consider mesons with different flavor, a frit-meson that we refer to, in analogy with it's four dimensional counterpart, 
as the "i3-meson," and a ud meson, the "pion." Flavors b, d, u will be denoted by a flavor index a = 1, 2, 3, respectively. 
While we refer to b as the "heavy" quark and to u and d as the "light" quarks, we do not make assumptions on the 
relative sizes of the masses, except when we consider the chiral limit for which we take m-i — — > (m u = md — > 0) 
holding mi (rrib) fixed. The form factors f± are defined bv|47| 

(n(p')\d(0)lA0)\B(p)) = (p + p r ),f + (q 2 ) (p-p'W-(<? 2 )- (41) 
where q = p — p' . Crossing this gives 

(0|d(0) 7| ,6(0)|S(p)7f(p')) = (p-p')Mq 2 ) - (p + p')J-(q 2 ). (42) 

where now q = p + p' . 

Explicit expressions for the form factors were determined in |13| by studying the crossed channel (|42|l . There it was 
shown that 

f+{q2) = ^ c^(y) (43) 

even n '* 

/ (q 2 ) = V Cn97rBn ^ q2 " > - V c ^Bn(q 2 )(ft 2 B -a4) / 44 n 

odd n n even n v 

where the label n refers to the tower of states with quantum numbers 12 = bd. The quantities g^Bm which are 
interpreted as triple meson couplings, are given by triple overlaps of wave- functions. In terms of Co — p_/q_, which is 
a function of q 2 only, it was found that, for even n 



g-KBn(q ) = T~T9n (45) 

while for odd n 

97rBn(q 2 ) = (46) 



where 



' r dz ^) { l)^ 2 ) { Z ^) Mz ) -If d z$( 13 )(4)^(f4)^(z) . (47) 

LO 1 — UJ LO ] r, LO 1 — LO 



jib - 
1 — LO 

As shown in Ref. j^l the form factors are super-convergent: at large q 2 



/± ~ | 9 2|l+/3 2 - ( 48 ) 
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By considering the integral § c dzf±(z) / '(z — q 2 ) over a contour C consisting of a circle at infinity deformed on the 
real line to avoid all the poles, one can show that the q 2 dependent couplings may be replaced by their value on-shell, 
e.g., 

Ml ) - 2. q 2 _ 2 - 2. ? 2_ M 2 (49) 
even n ,fc even n 

This is useful because one can show 13] that in the chiral limit, m 2 = m 3 — > 0+, one has g^s^M 2 ) ~ * f° r 71 7^ 0, 
and g^BB^s) ~ * ~ Mb/ctt where c,r = Cq 23 ^ is the normalized "pion" decay constant. Therefore, in the chiral limit 



fs/U 
1 - q 2 /^ 2 



Mf) = -fAf (so) 



B. Soft Functions 



Next we turn our attention to the non-diagonal analog of the shape function, the "soft function" defined, in analogy 
with the shape function, by 

F^J ^e-^-inip'Mx^Pie'So^b^Bip)) . (51) 

The soft function arises naturally in computations of B-meson decay amplitudes in SCET in which a pion is produced 
and is soft in the rest frame of the decaying J3-meson [l9l l33l l34j . 

The computation of the soft functions proceeds in complete analogy to that of the previous section. The first 
contribution to (|51[l. shown in Fig.0] is similar to what we found for the shape function above in the region < fc_ < 
p'_. But now there is an important difference, the first graph in Fig. 0] also gives a non- vanishing contribution from 
the region p'_ — P- < fc_ < 0. This new term cancels against a similar term from the second graph in Fig0] The 
result is 

F = ^ (1 + £ - ^ (1 - i/ U )9{ U - 00(0 + E ^ 12)(1 + ^ 2 )g " (g2) g(l + £ - (52) 

We have introduced lu = p'_/p^ and assumed to < 1. The meaning of £ is clear, £ = The functions (jj 1 ^ 13 ^ and 

0(32) correspond to the incoming and outgoing B and 7r mesons, respectively. The masses \i n are for the tower of 
states composed of b and u quarks, and g n (q 2 ) are functions of q 2 — {p — p') 2 given by 

9n(q 2 ) = - /W 13) (l - v)^ 2 \v/co)^ 2 \\^) - - J- f\v^\l - v)^\v/^\\^) (53) 
w Ju 1 — w J. — w J 1 — w 

It is straightforward to check that this result gives correctly the current form factor. Upon integration over £: 

(ttCpOI^V^ 15 \B(p)) = P-^ 2 ) + (p - p')- E £! 2 ^? ( 54 ) 

„ q M ™ 

where g„ was defined in l|53|) and k is given by 

«(? 2 ) = d^ (13) (1 + £ - ^)0 (32) (1 - , (55) 
Jo 

in agreement with the classical result of |38| ■ 

It will be useful to consider the crossed channel soft function. This is interesting in it's own right, but more 
importantly, it will be needed to investigate the chiral limit of the soft functions. The graphs are again given by 
Fig. 01 except with the direction of the pion momentum reversed. By direct computation we find 

F-=-J2 M f U f ] 9(x)9(l - x) , (56) 

1 Mn 



10 



where x is the momentum fraction of the outgoing light quark, x — 1 + cD£ = — +p'_), and g n ((! 2 ) is given by 

Eq. (|47|l . It is straightforward to check that this is in fact the analytic continuation of the last term on the right hand 
side of Eq. Q52JI . that is, the soft function in the region — (p_ — p'_) < fc_ < 0. 

We pause to make an interesting, though peripheral observation. The soft function in Eq. I|56[) vanishes for k- > 0, 
while the function in Eq. (|52|1 is explicitly non-vanishing in that region. Clearly one is not the analytic continuation 
(over q 2 for fixed k-) of the other. This violation of crossing symmetry is not particularly worrisome. Crossing 
symmetry is not a fundamental property of quantum field theories (and is not a "symmetry" ) . It is established on 
a case by case basis (see, e.g., |40j,|4lJ). It is surprising to see it fail because one has learned from experience that 
without exception one can show it. Except in the case at present. We suspect this case does not satisfy crossing 
symmetry because we are not considering an S-matrix element or a matrix element of a local operator but rather 
the matrix element of a non-local operator partially Fourier transformed. As soon as one integrates over fc_ crossing 
is recovered: one obtains the form factors l|44l) and (|54|l that can be easily checked to be the analytic continuation 
of each other. Alternatively one may simply state that, in this case crossing symmetry is the agreement between 
the analytic continuation of this two expressions for fc_ < 0, and that crossing symmetry has nothing to say about 
fc_ > 0. 

The computation of F + is straightforward. Since "07+^ = V'lV'z,! use °f Eq. Q indicates that the result is the 
same as in (|56H but with an additional factor of raim2/4pi_p2-, where pi.2 are the momenta of the quarks in Fig.0] 
Hence we obtain, in the "crossed" channel, 

2x(l - x) ^ q 2 - (il 

where x is the momentum fraction of the outgoing light quark, x = l+il>£; = — fc_/ (p_ +p'_). It is now easy to check 
that this give s the current form factor / + upon integration over x. One only needs to use Eq. i|31|l and the result, 
from Ref. [13J, that 

V = . (58) 



C. Chiral Limit 



Since the form factors simplify tremendously in the chiral limit, it is natural to ask whether this is also the case 
of the soft functions. Moreover, the chiral limit is well understood in four dimensions, and the methods used in 
understanding the chiral limit in four dimensions should apply as well in two dimensions. 

Consider the function 

„ r- K 

which appears in 1)56(1. We show in the Appendix that for fixed x this sum is super-convergent, that is 

G ( X ' ~ ( g 2)(l+fa) ' aS I 2 ~* 00 ■ ( 60 ) 

Hence one may integrate this in the complex q 2 plane over a closed contour that avoids the positive real axis and 
closes on a circle at infinity to show, just as for the form factor, that the numerators can be replaced by residues: 

G{x,?)=Y i ^f 9n ^\ (61) 

„ T- Mn 

We can now use the result that in the chiral limit the couplings to higher resonances vanish to simplify the soft 
functions of Eqs. (|56|l and l|57|l . as follows: 



(62) 
(63) 
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We have retained the vanishingly small mass ttt-2 in the expression for F + above. The chiral limit should be understood 
as having arbitrarily small mass, but care must be exercised in taking the limit m 2 = 0. For example, 777.2 plays a role 
regulating the behavior of F + as x — > 0. One can integrate over x to obtain a form factor, and then take the limit 
TO2 — * smoothly. 

To make this result completely explicit it is necessary to understand the chiral limit of 50 (Mb)- There is a subtlety 

here that must be addressed with care. For small but nonzero "pion" mass, \i v = ^q 2S \ the pole at q 2 = (1q falls 
bellow threshold. Therefore the value of the variable u) at q 2 = n 2 B is complex, and satisfies 

(l-a>) 2 = -4 + o(— ) (64) 
fJ-B VMs/ 

Comparing the exact expression for the form factor / + , Eq. I|49l) . with the result in the chiral limit, Eq. 1)51)0. and 
using the explicit form of the triple boson coupling given in Eqs. (|45|l - (|47l) . we obtain the chiral limit for go: 

3 ( M |)^-^I(1-^). (65) 

One can then check consistency by examining the chiral limit of the second form factor, /_. A similar computation 
yields 

M) - ,f* -x , (66) 

which is seen to coincide with Eq. H65[l when use is made of the resonance value of lj as given in Eq. (|(>4|) . 

Collecting results and expressing them in terms of the light-cone wave-functions, we see that, in the crossed channel, 



P - ^ , (67) 

U p- + pi q 2 - n B 



F + ^-2,l^^±4, (68) 

Jtt p+ q 2 - n B 

Finally, we can use these results to find the chiral limit of the soft functions in the "normal" channel. Since F_ in 
the crossed channel as given by Eq. I|56fl is the analytic continuation of the second term in Eq. I|52|) , we may express 
the soft function in the "normal" channel simply as 

2;U 2/g P- te-jy) _ (69) 
3 U P-- P- {P - P 1 ) 2 - Ms 

for — (p_— p'_) < k_ < 0. Here y is the argument of </4 12 ^ in (JS^Jl, namely, y = 1+^/(1— w) = (p_— pl+fc_)/(p_— p'_), 
which is interpreted as the momentum fraction carried by the b quark. Moreover, for < fc_ < p'_ we had the result 

0( 13 )(l + ^- w )^ 32 )(l-e/w) (70) 

which can be slightly simplified in the chiral limit, since the pion wave- function simplifies in that limit, (f)^ — > 1. 
Combining partial results we have as our final result the chiral limit of the soft function in the normal channel: 



2 P- <t>B-{y) 



P- -P- IP-P ) - MS 



(71) 



Here, we remind the reader, u = p'_/p_, £ = k-/p-, y = 1 + £/(l — £*;) = (p_ — p'_ + fc_)/(p_ — are the relevant 
momentum fractions and we assumed u> < 1. 

VI. DISCUSSION 
A. Shape Functions 

The result for the shape function, Eq. i|38fl , is rather surprising. The shape function "factorizes" into the product of 
light-cone wave functions for the incoming and outgoing mesons, but there is no reason to suspect this factorization 
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a priori. Certainly naive factorization would not have given this result. Naive factorization prescribes that one is to 
insert the vacuum in all possible ways into the operator that is being sandwiched with meson states, and extract the 
color singlet component only. But in this case the operator is a non-local quark bilinear. If one were to insert the 
vacuum between quark fields the resulting operator would have to be set to zero by this prescription, since it has no 
color singlet. Worse, vanishing of this quantity is protected also by conservation of quark number. 

There is little guidance for phenomenological models of the shape function. If one could only proof that the 
factorization observed in this paper is a more general consequence of the large N c limit, beyond the simple 1+1 
dimensional case, it would be interesting to use this as a first approximation to the shape function, with corrections that 
vanish in the large N c limit. Alternatively, since the shape function is directly measurable (up to power corrections) 
one could use this relation to infer the light-cone wave-function, which appears in many calculations but cannot be 
measured directly Brave souls may pursue this approximation as a conjecture, awaiting a proof in 4 dimensions: 

(72) 



/B-(fc-) = K \i> B -(A-k-)\ 0(A-fc_) 

where /s and ipB stand for the shape and light-cone wave-functions, A = Ao = (is — mi, and k~ 

This relation may give interesting constraints on both functions since both are separately constrained both by 
theory (e.g., through moments) and phenomenologically. For example, consider popular models for the light-cone 
wave-function and shape function of Refs. and |43j|. respectively: 



n &-(*-)■ 



(model) 
B— 



~K 



exp 



/■(model) 
JB- 



32 



(A - k^) 2 exp 



4 



k- 
~K 

(A-k.f 



9 (A 



(73) 
(74) 



where K and A are positive constants. While clearly these functions fail to satisfy the relation t[7'2) . the modifications 
needed to bring them into compliance with l|72|l are mild (since both are exponentials times polynomials). Whether 
this can be done while still satisfying general guiding principles in Refs. an d US used to propose these models is 
a question we do not embark on here. 



B. Soft function 

It was shown in [l9l | that the soft functions can be expressed in terms of the light-cone wave-functions in the chiral 
limit. Let us review the results from that paper. 

We begin by reviewing the effective Lagrangian for heavy mesons interacting with low energy pseudo-Goldstone 
bosons, the so called Heavy Meson Chiral Perturbation Theory (HMCHPT) EH ■ The effective Lagrangian, 
adapted to the two dimensional theory, is written in terms of a heavy meson super-field, 

H ° = - (-^) ^B a , (75) 





( 




7T+ 


K+ 


M = 




TT 




K° 






K~ 


K° 





2 

where B a = (&„, bd, B s ) are the positive energy fields for the B mesons, and a matrix of pseudo-Goldstone bosons, 

(76) 

We are considering the case of three light flavors, since this is what is phenomenologically relevant. Our computations 
in previous sections apply since the third light flavor does not contribute in any of the processes we have considered 
(to leading order in 1/N C ). The matrix of pseudo-Goldstone bosons appears in the Lagrangian through £ = e lM ^ 
(and S = £ 2 ). Under SU(3)l X SU(3)r these transform as, 

S -> LSi? t , £ -> L£[/t = U£R f and H a -> H b , (77) 

where the transformation U, defined through Eq. (|77|l . depends non- linearly on the pseudo-Goldstone bosons. In 
terms of these fields, the effective Lagrangian is, to lowest order in a momentum expansion and 1/m.b, 



f2 

C = -2i Tr [H {Q]a Vll d^H a Q) ] + Tr (d^Zd^) + A Tr [m,E + m q ^] 



i Tr [HM°HMv» + id^)J -f Tr [H a H b7vl5 - tP?)J +■■■ (78) 




FIG. 5: Non-resonant (left) and resonant (right) contribution to the soft function in Heavy Meson Chiral Perturbation Theory. 
The circled crosses denote insertions of the operators Ol,r. of Eq. 1841 . 

The coupling constant g determines the coupling of B mesons to pseudo-Goldstone bosons. The n—B — B interaction 
term in the Lagrangian is 



f ba D a 



(79) 



This is to be compared with the coupling g^BB which corresponds to a Lagrangian interaction term between three 
pseudo-scalar mesons 



gijke fl ud^Vid"(p j ipk 



(80) 



Here we have used the notation of Ref. 13] . In particular, we have a dimensionfull normalization constant relating 
the coupling with a carat and the couplings in Eqs. (|45|l and 1460 . g^BB — itib9kBB- Setting <fi(x) — tt(x), <fj[x) = 
e~ mihV ' x B(x) /y/m,b and (fk(x) — e imBV ' x B* (x) / y/mt, we see this is exactly the same interaction as in the chiral 
Lagrangian with the identification 



JttBB = m B g„ B B = m B 



f 



(81) 



Below we use this relation together we the 't Hooft model result c^g^s b = — 1 to eliminate g from the final expression 
for the soft function; see Eq. l|5rj|) . 

Operators in the theory are also constrained by the symmetries. For example, the left handed current L v a = q a lvPhb 
can be expressed in the low energy effective theory as |44ll45ll46j 



LI = -iaTrWP L H b tl} 



(82) 



The proportionality constant a can be fixed by requiring the correct matrix element between a i3-meson state and the 
vacuum, a — y/msfB (we have adopted here non-relativistic normalization of the B states, as is standard practice in 
HQET). 

Following Ref. we a PPly this technique to the non-local operators in the definition of the soft functions. To 
follow the notation in more closely, we introduce the null vectors n M = (1/V2, 1/V2) and — (l/y/2, —l/y/2). 



These have rv 



1 and n 



0. The non-local operators are 
' da; 



-ix k- 



-ix k- 



T 



q a {x-ytf A p Ll »b(Q) 



(83) 



which transform under the chiral group as (3^,1) and (1,3/?), respectively. The corresponding operators in the 
effective theory then are of the form 



0£ a (fe_) = Tr[a R (k.)P L ^HUba] 



(84) 



Using the fact that fiH — H, the most general form of ul.r is a linear combination of i/b and rfi. Matching to the 
matrix element between the vacuum and a B meson we obtain 



= a R (k-) = -ifB/\/mB[4^B- + $4>B+] 



(85) 
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It is now a straightforward exercise to compute the soft functions in the effective theory. To compare with our 
earlier results for F_ we consider the sum of operators Ol + O r . There are two contributions, as shown in Fig.JSJ the 
first is a contact term obtained by retaining a pion field in the expansion of £f, a and the second is a pole term from 
the no pion term in £& a . We obtain 



iy/mE(w(p')\(0 L + R )-\B(p)) = !f-^ B _(-k_) 



1 



, ty^u p 

' v ■ {p-p') 



(86) 



This is to be compared with our result, Eq. (|71|l . For the comparison we substitute m^v + p for p, so the meaning of 
p becomes the residual momentum, and retain only the leading term in an expansion in inverse powers of the heavy 
mass: 



V- 



(87) 



Using V- = u+ = 1/V2 and neglecting p' + = 2(j%/p'_ we see that except for the argument of the light-cone wave- 
function and the presence of 9- functions the two terms in (|86[1 agree with the corresponding terms in (I87|l . 

The result is somewhat puzzling, since the form of Eq. i|86|) was obtained through general arguments in Ref. |19| . 
However, the arguments presented there are not quite rigorous. We can partially resolve the discrepancy as follows. 
Consider gauging the vector flavor symmetries of the QCD Lagrangian. We gauge both the full, non-abelian, flavor 
symmetry of the light quarks and the flavor symmetry of the heavy quark, "6-number." The gauge fields are taken as 
background fields and are introduced solely for the purpose of maintaining local invariance, but will be set to zero at 
the end of the calculation. Now, instead of the integrated operators in l|83|) . we consider the un- integrated, non-local 
operators 

o£ a ( x -,o) = T(q a (x~y^ A p B Yb(o), 

Under a local vector flavor symmetry transformation these transform as 

Ol R ^B(0)O»V(x-y, 



(89) 



where B and V stand for the unitary transformations of 6-number and light flavor, respectively. The corresponding 
operators representing these in the chiral Lagrangian should transform the same way. Hence we write 



O Ra (x-,0) =Tr[a R (0,x-)P L ^H c (0)P cb (0,x~)€ b a(x~)]- 



(90) 



Here otL,R are coefficient functions, to be determined, and /3(0,x ) = Texp(— i L A v ), where is the background 
vector gauge field associated with the light quark flavor vector symmetry, transforms as 



(3(0, x~) -> V(0)(3(0,x~)V\x- 



(91) 



Turning off the background fields, Fourier transforming these operators and repeating the argument above we now 
obtain 



t^(n(p')\(0 L + O r )-\B(j>)) = t± 



$B-(pL - k^)6(p'_ - fc_) + 2&,_(-A;_)0(-fc_) 



eu^p' v 



v-{p- p') 



(92) 



This only fails to reproduce Eq. (|87|l in that the first term (the "contact" interaction) is missing the (9-function enforcing 
fc_ > 0. In the 't Hooft model calculation the origin of this restriction is clear: the contact terms involves the product 
of the wave-functions of the B- and 7r-mesons, and the latter requires fc_ > 0. However, the chiral Lagrangian operator 
we have constructed has lost this information. It would be interesting to see how this problem could be fixed. It seems 
likely that some understanding of the representation of non-local purely light-quark operators in the chiral/SCET 
Lagrangian is necessary, since this should carry information about the pion wave-function. However, since the physical 
basis for the restriction that fc_ > is pretty clear, it does make sense to adopt this result in the full 4-dimensional 
analysis, now modified by replacing the non-local operators (|9t)|l for l)84[l. 
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VII. CONCLUSIONS 



We have computed the analogues of the shape function and soft functions of 4-dimensional HQET/SCET in QCD 
in 1+1 dimensions (the '"t Hooft Model"). These are matrix elements between mesons of non-local operators. 

Our main results are as follows. First, the shape function, /s(fc_), is, up to a multiplicative constant fixed by 
normalization conditions, the square of the the light-cone wave function, 

/ B (fc_) = 2 e(-k-) . 

Second, in the chiral and heavy meson combined limits the soft function for B to it transitions (defined in (151(1 ) is 
completely determined by the B meson light-cone wave-function, as argued on general grounds in Ref. [ijfl . However, 
we find that the representation of the non-local operator in the effective theory presented in Ref. is not correct. 
A better representation of this operator is given by a non-local operator in the effective theory, e.g., 

Ol a (x-,0) = Tr[a L (O,x-)P R rH c (O)0 cb (O,x-)^ a (x-)] . 

This operator almost correctly reproduces the result from direct computation in the 't Hooft Model. It only fails 
to reproduce a 8- function restricting the momentum of the light quark (see Sect. IVTBl for details), which however is 
easily understood on physical grounds and can be adopted (albeit in a somewhat ad-hoc fashion) in computations in 
four dimensional QCD. 

We cannot argue that the first of these results is applicable in four dimensional QCD. However, the two results are 
tightly connected in the 't Hooft model. This suggests that there may in fact exist a way to justify the first result 
in four dimensions. It is interesting to note in this regard that the same argument that led us to consider non-local 
operators in the effective theory for the soft-functions, gives that the effective operator in the effective theory for the 
shape function is also non-local, roughly 

Or ~ Tr [a(x-,0)H a (x-)(3(x-,0) ab TH b (0)] . (93) 

Hence, if one could argue that a(x~,0) factorizes, then the result / cx ip 2 would follow automatically. Were this 
established, interesting phenomenological constraints relating exclusive to inclusive £?-decays would follow. Clearly 
this is a topic that deserves much further investigation. 
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APPENDIX: SCALING THEOREMS 

We follow the work of Einhorn closely 38]. Our scaling functions are similar but not the same as his. Let 

^)(x,, ;g 2 )=E^ ! ^V M - (A.1) 

q 

One then shows this has a limit, which defines the scaling functions 

h£\t;y) = lim q 2 g^(±\,y;q 2 ) (A.2) 

q^^oo q 

Note that integrating h(£;y) over y gives Einhorn's scaling function h(£). The scaling function satisfies 

l-f)4 1) (^) + /d, 7 |^ = %) ! (A.3) 
from which it immediately follows that 

h£>(Z;y)~^ as £ -> 0, fixed y. (A.4) 
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Knowing that Q has a finite limit one can determine the large q 2 scaling of G(x, q 2 ) in Eq. (|60(l . To this end write 
explicitly the first term in g n (q 2 ) from Eq. H47|) in G(x,q 2 ) as 



i-cj dz ^ ( o J$ 4-o j ^ ? 2 -^ 



q r\ r\ ^{v^X^g^iuv.x^ 2 ) 

' av at 



1-0 



o Jo 



t + 



T^(l--))' 



1 ,1/0 0(13) ( w )^(32) (0^1^(12)^ g..^ 



dw / d«- — V^H^ — (A.5) 

Ji {u-vy 

In the second line we have changed integration variables (v = z/0), used the integral representation of and expressed 
the sum in terms of the function Q of IjA.lfl . The third line uses a further change of variables uju = t(l — 0) + uj so 
the denominator factors, 

((1 - tb)t + 0(1 - v)) 2 = uj 2 (u - v) 2 , (A.6) 

which leaves the integrand in a form suitable for exploration of the small behavior. This corresponds to large q 2 , 
q 2 ~ m%/u). 

Using 0^ 32 ) (x) w cx@ 3 as x — > 0, the behavior of (|A.5(1 is obtained by writing it in terms of the scaling function: 

c (m|)A f 1 ^ n 4M(v)(u-l)P»h¥\mlv;x) ?) 



Jo~"Ji~~ (»-'') 
Analogous steps are taken for the second term in G(x; g 2 ) 

, z - O ^ 0l 12) (a;)0l 12) (z) 



1 /L$( 13 )(4)0 (32) (^)V- 

wis y l-u>^ 



q 2 - Mn 



1-0 f 1 ^ f dt ^(t)^(v)g^(Q + (l-ij)v,x-, q 2 ) 



O 7o Jo (t-l-l^£ w ) 2 
= — at au 



Jo Si 



{u - t) 2 

4 ^ ^ 13 )(t)(«-l)ft^ 12) (m| U ;x) _ 
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